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ABSTRACT 

The presence of hot gaseous coronae around present-day massive spiral galaxies is a fundamental 
prediction of galaxy formation models. However, our observational knowledge remains scarce, since 
to date only four gaseous coronae were detected around spirals with massive stellar bodies (> 2 x 
10^^ Mq). To explore the hot coronae around lower mass spiral galaxies, we utilized Chandra X- 
ray observations of a sample of eight normal spiral galaxies with stellar masses of (0.7 — 2.0) x 
10^^ M0. Although statistically significant diffuse X-ray emission is not detected beyond the optical 
radii (~20 kpc) of the galaxies, we derive Scr limits on the characteristics of the coronae. These limits, 
complemented with previous detections of NGC 1961 and NGC 6753, are used to probe the Illustris 
Simulation. The observed Scr upper limits on the X-ray luminosities and gas masses exceed or are at 
the upper end of the model predictions. Eor NGC 1961 and NGC 6753 the observed gas temperatures, 
metal abundances, and electron density profiles broadly agree with those predicted by Illustris. These 
results hint that the physics modules of Illustris are broadly consistent with the observed properties 
of hot coronae around spiral galaxies. However, a shortcoming of Illustris is that massive black holes, 
mostly residing in giant ellipticals, give rise to powerful radio-mode AGN feedback, which results in 
under luminous coronae for ellipticals. 

Subject headings: cosmology: theory - galaxies: spiral - galaxies: ISM - X-rays: galaxies - X-rays: 
general - X-rays: ISM 


1. INTRODUCTION 

The existence of gaseous X-ray coronae in the 
dark matter halos of massive gal axies is predic t ed by 
galaxy formation models (e.g. iWhite fc Reeij Il978t 
White fc Frenklll99lHToft et al. I2002t iRasmussen et akl 

2009t iCrain et al.l l201Cl[ iVogelsberger et akl 120131 : 

Marinacci et al. IMi. Although simulations agree that 
these hot coronae should be ubiquitous, various models 
made very different predictions. The disagreement 
among models originates from two sources. First, the 
incomplete realization of the physical processes may 
result in different galaxy populations. Specifically, many 
earlier models lacked efficient supernova and/or AGN 
feedback, which often yielded overly massive stellar 
components relative to their dark matter halos due to 
the strong gas cooling. Second, using different hydrody- 
namical schemes resulted in large sy stematic differences 
in th e global state of baryons (|Vogelsberger et akl 
1201211 . Thus, the predicted properties of hot coronae 
also sen sitively depend on th e applied hydrodynamical 
scheme (iBogdan et al.l[^13ajl . The ever increasing com¬ 
puting power plays a major role to overcome the issues 
related to incomplete modeling of the physical processes. 
Indeed, state-of-the-art simulations are now capable 
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of studying large cosmological volumes combined with 
high numerical resolution, which simultaneously allows 
the study of large galaxy samples a nd robust modeling 
of com plex physical processes (e.g. IVogelsberger et aH 
l2014ajl . Additionally, detailed tests helped to under¬ 
stand the differences between the recently developed 
moving mesh and the traditional smoothed particle 
hydrodynamics codes. 

Although hot X-ray coronae are of fundamental im¬ 
portance, their detection and characterization poses ma¬ 
jor observational challenges around spiral galaxies, which 
system offer a clean test of galaxy formation models 
((Bogdan et al.1 [^13all . The main difficulty is due to 
the faint but extended nature of the coronae around spi¬ 
rals, which in turn results in low signal-to-noise ratios. 
Additionally, the brightest X-ray coronae are associated 
with the most massive galaxies, which systems are rare 
in the local Universe {D < 100 Mpc). Despite these 
challenges, numerous observational studies aimed to ex¬ 
plore hot gaseous coronae around nearby spiral galax¬ 
ies, which system s offer a clean test of galaxy forma¬ 
tion models (e.g. iBenson et al.l 120001 : IRasmussen et al.1 
l2009f l. However, the first X-ray coronae around spirals 
have only recently been detected using Chandra, XMM- 
Newton, and ROSAT data, which revealed gaseous 
coronae around NGC 1961, UGC 12591, NGC 6753, 
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TABLE 1 

The properties of the sample galaxies. 


Name 

Distance 

(Mpc) 

( 1 ) 

1 ' scale 
(kpc) 
( 2 ) 

Nh 

(cm~^) 

(3) 

Lk 

(Lk,©) 

(4) 

M*/iK 

(Mq/Lk,©) 

(5) 

M* 

(Me) 

( 6 ) 

SFR 

(Me yr-i) 
(7) 

Morph. 

type 

( 8 ) 

r 200 

(kpc) 

(9) 

F>26 

(’) 

( 10 ) 

NGC 266 

60.3 

17.55 

5.7 x 10^0 

2.5 X 10“ 

0.80 

2.0 X 10“ 

2.4 

Sab 

410 

2.95 

NGC 1097 

14.2 

4.12 

2.0 X 10^0 

1.3 X 10“ 

0.77 

1.0 X 10“ 

5.8 

SBb 

238 

10.47 

NGC 2841 

14.1 

4.09 

1.3 X 10^0 

1.5 X 10“ 

0.78 

1.2 X 10“ 

0.8 

Sb 

363 

6.92 

NGC 5005 

18.5 

5.38 

1.1 X lO^o 

1.9 X 10“ 

0.77 

1.4 X 10“ 

4.9 

SABb 

273 

4.79 

NGC 5170 

21.6 

6.29 

6.7 X lO^o 

8.5 X 10i° 

0.83 

7.1 X lOio 

0.4 

Sc 

269 

7.94 

NGC 5529 

36.0 

10.45 

1.0 X 10^0 

1.0 X 10“ 

0.75 

7.8 X 10^° 

2.0 

Sc 

314 

5.75 

NGC 5746 

19.8 

5.79 

3.4 X 10^0 

1.4 X 10“ 

0.77 

1.1 X 10“ 

0.6 

SABb 

353 

7.24 

ESO 445-081 

69.3 

20.16 

5.0 X 10^0 

1.2 X 10“ 

0.78 

9.4 X 10i° 

5.6 

SBbc 

590 

1.70 


Note. Columns are as follows. (1) Di stance taken fromfrullv et aO 1)20091) (2) 1' scale at the adopted distance. (3) Line-of-sight column 
density provide d by the LAB su rvey JKalberla et al.ll2005l) . f4) Total K-band luminosity from 2MASS. (5) K-band mass-to-light ratios 
computed from IBell et al.l <120031) using the B — V color indices of galaxies. ( 6 ) Total stellar mass based on the K-band luminosity and the 
K-band mass-to-light ratios. (7) Star formation rate computed from the IRAS 60 fim and 100 fim flux - for details see Section 13.21 ( 8 ) 
Morphologica l type, taken from HyperLeda. (9) Virial radius of the galaxies, estimated from the maximum rotation velocity following 
IBogdan et al.l (I2013al) . (10) Major axis diameter of the D 25 ellipse. 


and NGC 266 (I Anderson et alJ 120111 : iDai et all 120111 : 
IBogdan et alJ I2013al lbfl. Due to the large distances of 
these galaxies and/or the relatively shallow X-ray ob¬ 
servations, the detailed characterization of the coro- 
nae was only p ossible for NGC 1961 and NGC 6753 
(|Anderson et al.l[2011l: IBogdan et al.ll2013all . 

Despite these advances, the presently available sam¬ 
ple is still scarce and only includes massive galaxies 
(M* > 2 X 10^^ Mq), which presents limitations in using 
the hot coronae as probes of galaxy formation models. 
While extra-planar diffuse emission has been detected 
around several lower mass spiral galaxies (|Li et al.ll2006L 
120111: IBogdan fc GilfanovI I2008D . this gaseous emission 
cannot be attributed to the hot corona of infalling gas. 
Indeed, the diffuse gas in these galaxies do not extend 
to large radii, and exhibits a bipolar morphology. The 
morphology combined with the mass and energy budget 
of the galaxies hint that the gas is outflowing and the 
outflows are sustained by the energy input from type la 
supernovae. Additionally, several studies have explored 
the presence of diffuse emission within the central re - 
gions of spiral galaxies (jLi fc Wangll2013HLi et al.ll2014l l. 
Although a certain fraction of this emission may origi¬ 
nate from the infalling hot gas, it is nearly impossible to 
disentangle this emission from other diffuse X-ray emit¬ 
ting components, in particular from the hot gas associ¬ 
ated with star-formation. Therefore, it is desirable to 
probe hot coronae beyond the optical radii of the galax¬ 
ies. Exploring extended X-ray coronae around lower 
mass spirals will allow us to place stronger constraints 
on galaxy formation simulations and to probe a broader 
galaxy population. However, due to their faint nature, 
the X-ray coronae may remain undetected around galax¬ 
ies with AL < 2 X 10^^ Mq. While non-detections only 
allow limits on various properties of the coronae, these 
limits provide invaluable input to probe galaxy forma¬ 
tion models. Thus, the goal of this work is to extend the 
presently available sample with a focus on lower mass 
spiral galaxies. 

The Illustris Sim ulation, which is base d on the moving- 
mesh code AREPO (iSoringel et al.l[20lRl . offers the ideal 
framework to compare the observational results with 


those predicted by a modern hydrodynamical cosmolog¬ 
ical simulation. Illustris includes the main physical pro¬ 
cesses that influence galaxy evolution, such as primor¬ 
dial and metal line cooling, star-formation, chemical en- 
richment, energetic feedbac k from supernovae and AGN 
(|Vogelsberger et aT]l2014bll . Therefore, Illustris results 
in a real istic galaxy population o f elliptical and spiral 
galaxies (jVogelsberger et all 120131 1. Additionally, Illus¬ 
tris was performed in a large cosmological box, which 
consists of a broad population of spiral galaxies with 
M* > 10^^ Mq, allowing a comparison with observa¬ 
tional results. Thus, the goal of our study is to probe 
hot X-ray coronae around spiral galaxies and compare 
the observed properties with those predicted by Illustris. 

The paper is structured as follows. In Section 2 we 
introduce the observed and simulated galaxy sample. In 
Section 3 the data reduction is described. We present the 
results in Section 4, where we discuss the constraints on 
the characteristics of the hot X-ray coronae. In Section 
5 we compare the observed and predicted properties of 
the coronae. We discuss our results in Section 6 and 
summarize in Section 7. 

2. GALAXY SAMPLES 
2 .1. The observed galaxy sample 

To probe the hot X-ray coronae around spiral galaxies, 
we utilize Chandra X-ray observations. While XMM- 
Newton offers larger effective area, a particular advan¬ 
tage of Chandra is its lower and more stable instrumen¬ 
tal background level, which makes it ideal to study the 
faint X-ray emission expected from hot coronae around 
spiral galaxies with M* < 2 x 10^^ Mq. 

To this end, we searched for spiral galaxies with pub¬ 
licly available Chandra observations. To ascertain the 
morphological type of the galaxies, we relied on the Hy¬ 
perLeda data base and included Sa-Sc systems. We 
further filtered the sample by selecting galaxies in the 
distance range of 14 < D < 40 Mpc. The lower dis¬ 
tance boundary assures that at least the inner regions of 
the coronae can fit in the Chandra field-of-view, while 
the higher boundary was introduced to maximize the 
signal-to-noise ratio. We also required that the sample 
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Fig. 1.— DSS B-band images of the spiral galaxies in our sample. All galaxies reveal late-type morphologies and have a prominent disk 
component. The images show that the selected galaxies are relatively undisturbed, only NGC 1097 may be experiencing a minor merger. 
The sample galaxies exhibit a wide range of inclinations, including face-on and edge-on systems as well. 


galaxies be sufficiently massive (M* > 7 x 10^° Mq), 
which ensures that they reside in a massive dark mat¬ 
ter halo, and host a more luminous and hotter X-ray 
corona. Additionally, we excluded those galaxies, which 
are located in rich galaxy groups/clusters, and only re¬ 
tained field galaxies or systems residing in poor groups 
with < 20 members. This assures that we can probe the 
hot corona associated with the galaxy and not a larger- 
scale group/cluster atmosphere. Moreover, we excluded 
galaxies undergoing a significant merger event, since the 
dark matter halos, and hence the hot gaseous coronae, 
of these galaxies may be disturbed. Finally, we removed 
starburst galaxies (for details see Sectional) and systems 
with high star-formation rates (> 10 Mq yr“^), since 
these galaxies could be prone to starburst driven winds. 
Based on these criteria we identified six nearby spiral 
galaxies. 

This sample was further expanded by two additional 
spiral galaxies, NGC 266 and ESO 445-081. While 
these galaxies are more distant, their ROSAT observa- 
tions indicate the pr esence of luminous X-ray coronae 
(jBogdan et al.l2013b( l. Thus, we utilized follow-up Chan¬ 
dra observations to further study these galaxies. 

Figure [1] shows the Digitized Sky Survey (DSS) B-band 
images of the 8 spiral galaxies in our sample. While 
all galaxies reveal late-type morphologies, they exhibit a 
wide range of inclinations from face-on to edge-on sys¬ 
tems. The physical properties of these spirals are listed 
in Table [TJ 

2.2. The simulated galaxy sample 

To compare the observed characteristics of the coro¬ 
nae with those predicted by Illustris, we used two ap¬ 
proaches. First, we selected all galaxies independently 


TABLE 2 

The list of analyzed Chandra observations. 


Galaxy 

Obs ID 

Toha (ks) 

ffiilt (ks) 

Instrument 

NGC 266 

16013t 

85.0 

85.0 

ACIS-I 

NGC 266 

16301 

41.5 

41.5 

ACIS-I 

NGC 1097 

161lt 

5.4 

5.0 

ACIS-S 

NGC 1097 

2339 

5.7 

5.0 

ACIS-S 

NGC 2841 

389 

1.8 

1.8 

ACIS-S 

NGC 2841 

6096+ 

28.2 

21.3 

ACIS-S 

NGC 5005 

4021 

4.9 

4.7 

ACIS-S 

NGC 5170 

3928 

33.0 

33.0 

ACIS-I 

NGC 5529 

4163 

89.2 

71.0 

ACIS-I 

NGC 5529 

12255 

60.4 

50.8 

ACIS-S 

NGC 5529 

12256+ 

118.5 

96.0 

ACIS-S 

NGC 5529 

13118 

44.6 

38.3 

ACIS-S 

NGC 5529 

13119 

54.3 

48.0 

ACIS-S 

NGC 5746 

4021 

36.8 

33.6 

ACIS-I 

ESO 445-081 

16302 

39.2 

39.2 

ACIS-I 


t The coordinate system of these observations was used as refer¬ 
ence when merging observations. 


of their morphological types, whose stellar mass is in 
the range of (0.2 — 20) x 10^^ Mq. Due to their shal¬ 
low potential well, the dominant fraction of galaxies with 
M* < 2 X 10^° Mq cannot confine a significant amount 
of hot X-ray gas. This approach offers a statistically sig¬ 
nificant galaxy sample, allowing us to characterize the 
gaseous coronae from low stellar mass systems to the 
most massive galaxies. Additionally, this large galaxy 
sample removes the distorting effects of outliers, which 
ensures that we can explore the general trends in the 
simulated properties of the hot coronae. 
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Second, we analyz ed a sub-sample of spiral galaxies 
that were selected bv iVoeelsberger et al.l (|2014al) . These 
42 “textbook” spirals are sufficiently massive and are disk 
dominated systems. Therefore, the “textbook” spirals 
are the analogs of the observed galaxy sample, and hence 
offer an ideal basis to confront the Chandra observations 
with the Illustris Simulation. 

3. DATA REDUCTION 
3.1. Chandra X-ray observations 

The analysis of the Chandra data was performed with 
standard CIAO tools (version 4.6 and CALDB version 
4.6.2). For each galaxy we analyzed all available ACIS-S 
or ACIS-I observations. The total combined exposure 
time of the analyzed data is 648.5 ks. The analyzed 
Chandra observations are listed in Tabled] 

The first step of the data analysis was to reprocess 
all observations using the CHANDRA_repro task, which 
results in enhanced data quality and better calibrated 
observations. After this we identified and removed time 
intervals that were contaminated with background flares. 
Since the detectors are most sensitive to flares in the 
2.3 — 7.3 keV band, we used this energy range to identify 
high background periods. For each observation, we ex¬ 
cluded those time intervals that were 2a above the mean 
count rate level. Removing the flares resulted in a ~ 12% 
drop in the exposure times, and a total net exposure 
time of 574.2 ks. However, note that the ACIS-S detec¬ 
tors are more susceptible to background flares, hence a 
larger fraction of contaminated time intervals were ex¬ 
cluded from these data (Tabled). 

In general, bright point sources can add a significant 
contribution to the truly diffuse gaseous emission, hence 
it is necessary to exclude them when studying diffuse X- 
ray emitting components. To detect bright point sources, 
we ran t he wavdetect task using the parameters de¬ 
scribed inlBogdln^^^nSaQ^ ([lOO^). This procedure re¬ 
sults in sufficiently large source cells that include > 97% 
of the source counts. These source cells were used to 
mask out bright point sources from the further study of 
the diffuse emission. 

To account for vignetting effects and convert the ob¬ 
served counts to flux units, we produced exposure maps 
assuming a thermal plasma model (apec in XSPEC) . The 
temperature of the model was 0.2 keV and the metal 
abun dances were set to 0.1 Solar (|Grevesse fc SauvO 
Il998fl . The particular choice of spectrum was motivated 
by the results of our earlier observations (| Bogdan et al.l 
l2013all and the Illustris Simulation. Indeed, Illustris pre¬ 
dicts similar gas temperatures for a hot corona surround¬ 
ing a galaxy with M* = 10^^ Mq, which is close to the 
median stellar mass of our sample. Additionally, we 
observed sub-solar metal abundances 0.1 Solar) for 
NGC 1961 and NGC 6753, which is in fair agreement 
with that predicted by Illustris (Section [O])- 

Since we aim to study the faint diffuse X-ray emission, 
it is crucial to precisely account for the background emis¬ 
sion. We used two different approaches to subtract the 
background components. For the more distant galaxies 
in our sample, we used local fields to subtract the back¬ 
ground level. Using internal background fields offers a ro¬ 
bust way to subtract all background components without 
the need to renormalize the background level. We note 


that due to the large virial radii of the galaxies in our 
sample, the applied background regions may be within 
the extended coronae of the galaxies. However, at large 
radii (> 0.3r2oo) the gas is expected to have very low den¬ 
sity (see Figure 0), implying that the extended coronae 
remain well below the detection threshold. To confirm 
this picture, we extracted surface brightness profiles of 
the sample galaxies in the 0.3 — 2 keV energy range, and 
found that the vignetting corrected profiles do not show 
variations at these large radii. This hints that at radii 
> 0.3r2oo the contribution of gaseous coronae are negli¬ 
gible compared to the sky and instrumental background 
components. For the more nearby galaxies in our sample, 
we relied on the AGIS blank-sky background files to sub¬ 
tract the background components. To ensure that soft 
X-ray sky background of the sample galaxies is similar 
to that of the background maps, we checked the fluxes of 
ROSAT All-Sky Surve y R4-R5 band ima^ s at the po¬ 
sition of the galaxies ((Snowden et al.l[T997ll . We found 
that the galaxies in our sample have similar soft X-ray 
brightness to those used for constructing the AGIS back¬ 
ground files. Due to the stable spectrum of the instru¬ 
mental background, the Chandra blank-sky files can be 
tailored to a specific observation. To this end, we normal¬ 
ized the blank-sky background files using the observed 
count rate ratios in the 10 — 12 keV band. The accu¬ 
racy of these background subtraction methods were com¬ 
pared for those galaxies, where both procedures could be 
applied. We concluded that the background subtracted 
data in these galaxies are in good agreement with each 
other, and the accuracy of background subtraction is a 
few per cent. 

3.2. Near- and far-infrared data 

The near-infrared images offer a robust means to trace 
the stellar light. Therefore, we rely on the K-band images 
of the Two-Micron All Sky S urvey (2MASS) to d erive the 
stellar mass of the galaxies (|Jarrett et al.ll200^ . Using 
the K-band magnitudes and the distances of the galaxies 
we derive the absolute K-band magnitudes and luminosi¬ 
ties fol lowing the procedure described in iBogdan et al.l 
(|2013ajl . The K-band luminosities are converted to stel¬ 
lar mass using the mass-to-light ratios computed from 
the corresponding B — V color indices an d results of 
galaxy evolution modeling (iBell et al.ll2003ll . 

The star-formation rates of the sample galaxies are 
computed from the 60 /rm and 100 /rm far-infrared fluxes 
obtained by the Infrared Astr onomical Satellite {IRAS) 
and the relation established bv iKennicuttI (jl998D . To de¬ 
rive the total far-infrared luminosity of the sample galax¬ 
ies, and hence derive the star -formation rates, we f ollow 
the methodology described in iBogdan et al.l (l2013all . We 
conclude that all of the sample galaxies have low star- 
formation rates relative to their stellar mass, implying 
that none of the galaxies is a starburst system. The stel¬ 
lar masses and star-formation rates of the sample galaxies 
are given in Table [1] 

4. RESULTS 

4.1. Searching for hot X-ray coronae 

To identify hot coronae around the sample galaxies, 
we employed four different methods. First, we searched 
for diffuse X-ray emission by visually inspecting the back- 
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TABLE 3 

Constraints on the parameters of the diffuse hot gaseous coronae. 


Galaxy 

Radial range 

^0.5-2keV,abs 

^'bol 

Afgas 

rie 

^cool 

M 


^"200 

erg 

erg s~^ 

Mq 

cm“° 

Gyr 

M© yr-° 

NGC 266 

0.05-0.15 

< 1.9 X 10'*° 

< 8.9 X 10"*° 

< 1.9 X 10°° 

< 8.0 X 10-* 

> 2.9 

< 6.5 

NGC 1097 

0.05-0.15 

< 1.5 X 10°° 

< 6.1 X 10°° 

< 2.2 X 10° 

< 5.0 X 10-* 

> 4.6 

< 0.5 

NGC 2841 

0.05-0.15 

< 2.9 X 10°° 

< 1.1 X 10“*° 

< 5.8 X 10° 

< 3.6 X 10-* 

> 6.4 

< 0.9 

NGC 5005 

0.05-0.15 

< 2.3 X 10°° 

< 8.7 X 10°® 

< 3.3 X 10° 

< 4.9 X 10“* 

> 4.7 

< 0.7 

NGC 5170 

0.05-0.15 

< 3.0 X 10°° 

< 1.6 X 10'*° 

< 4.3 X 10° 

< 6.6 X 10-* 

> 3.5 

< 1.2 

NGC 5529 

0.05-0.15 

< 4.7 X 10°° 

< 1.8 X 10"*° 

< 5.8 X 10° 

< 5.6 X 10-* 

> 4.1 

< 1.4 

NGC 5746 

0.05-0.15 

< 3.0 X 10°° 

< 1.3 X 10'*° 

< 5.9 X 10° 

< 4.0 X 10“° 

> 5.7 

< 1.0 

ESO 445-081 

0.05-0.15 

< 2.9 X 10'*° 

< 1.4 X 10"*° 

< 4.1 X 10°° 

< 6.0 X 10-* 

> 3.8 

< 10.9 

STACKt 

0.05-0.15 

< 2.0 X 10°° 

< 7.3 X 10°® 

< 3.2 X 10° 

< 3.8 X 10“* 

> 6.0 

< 0.5 

NGC 266 

0.15-0.30 

< 4.0 X 10'*° 

< 1.9 X 10"** 

< 7.3 X 10°° 

< 4.4 X 10“* 

> 5.3 

< 13.8 

NGC 1097 

0.15-0.30 

< 2.9 X 10°° 

< 1.2 X 10'*° 

< 8.4 X 10° 

< 2.6 X 10-* 

> 8.9 

< 0.9 

NGC 2841 

0.15-0.30 

< 1.7 X 10°° 

< 6.6 X 10°° 

< 1.2 X 10°° 

< 1.0 X 10“* 

> 22.4 

< 0.5 

NGC 5005 

0.15-0.30 

< 5.0 X 10°° 

< 1.9 X 10'*° 

< 1.3 X 10°° 

< 2.7 X 10“* 

> 8.6 

< 1.5 

NGC 5170 

0.15-0.30 

< 4.3 X 10°° 

< 2.2 X 10"*° 

< 1.4 X 10°° 

< 3.0 X 10-* 

> 7.8 

< 1.8 

NGC 5529 

0.15-0.30 

< 8.0 X 10°° 

< 3.0 X 10'*° 

< 2.1 X 10°° 

< 2.7 X 10“* 

> 8.4 

< 2.4 

NGC 5746 

0.15-0.30 

< 2.1 X 10°° 

< 9.0 X 10°° 

< 1.3 X 10°° 

< 1.2 X 10-* 

> 18.3 

< 0.7 

ESO 445-081 

0.15-0.30 

< 6.4 X 10'*° 

< 3.1 X 10'*° 

< 1.7 X 10°° 

< 3.3 X 10-* 

> 6.8 

< 24.4 

STACKt 

0.15-0.30 

< 2.7 X 10°° 

< 9.8 X 10°° 

< 1.0 X 10°° 

< 1.6 X 10-* 

> 14.4 

< 0.7 


^ Based on the stacked data of a sub-sample of six galaxies: NGC 1097, NGC 2841, NGG 5005, NGC 5170, NGG 5529, and NGC 5746. 



Fig. 2.— Stacked 0.3 — 2 keV band X-ray surface brightness 
profile of the diffuse emission for the eight galaxies in our sample. 
To construct the profiles, we used circular annuli centered on the 
optical centroid of each galaxy. Vignetting correction is applied 
and all background components are subtracted. Over plotted is the 
2MASS K-band surface brightness profile, which traces the stellar 
light distribution. The K-band light is normalized to match the 
last X-ray detected bin at 17 kpc. Note that the X-ray surface 
brightness for the last three bins are consistent with 0. 


ground subtracted and vignetting corrected Chandra im¬ 
ages in a soft (0.3 — 1 keV), a broad (0.3 — 2 keV), and a 
hard band (1 — 2 keV). Probing different energy ranges 
is useful, as it offers the opportunity to detect an ionized 
gas component with relatively low (< 0.3 keV) or higher 
(> 0.3 keV) gas temperatures. While a diffuse X-ray 
glow is associated with some of the spiral galaxies in our 
sample, this emission is focused in the innermost regions. 
However, within these regions several additional X-ray 
emitting components are present, such as unresolved low- 


mass and high-mass X-ray binaries, cataclysmic vari¬ 
ables, active binaries, young stars, young stella r objects, 
and diffuse gas associated with star formation (iGilfanovI 
2004 Ifevniytsev et al.l [2008t iBogdan &: GilfanovI 120111 : 
Mineo et al. IMf Given the overall difficulty to assess 
the importance of each of the above listed components 
to the diffuse gaseous emission, we must explore the hot 
coronae beyond the optical radii of the galaxies. 

As a second approach, we probed the existence of hot 
coronae by constructing radial surface brightness profiles 
for each galaxy in the 0.3 — 1 keV, 0.3 — 2 keV, and 1 — 2 
keV bands. To construct the profiles we utilized the back¬ 
ground subtracted and vignetting corrected data, and ex¬ 
tracted circular annuli centered on the optical centroid 
of each galaxy. The X-ray profiles are compared with 
the stellar light distribution inferred from the 2MASS 
K-band images. Based on these profiles we did not de¬ 
tect an extended X-ray emitting component beyond the 
optical radii. 

Third, to increased the signal-to-noise ratio of the data, 
we co-added (i.e. stacked) the individual surface bright¬ 
ness distributions in the three studied energy ranges. 
In Fig. [21 we show the stacked 0.3 — 2 keV band sur¬ 
face brightness profile of the diffuse X-ray emission as 
well as the stacked K-band light distribution. While the 
X-ray and K-band surface brightness distributions are 
markedly different, the stacked profile does not reveal a 
statistically significant diffuse X-ray emitting component 
beyond ~ 20 kpc. In agreement with these, the soft and 
hard band profiles also did not reveal the presence of 
diffuse X-ray emission beyond the extent of the stellar 
light. 

Finally, we probed the existence of an extended X- 
ray emitting component in the (0.05 — 0.15)r2oo and the 
(0.15 —0.30)r2oo regions in the soft, broad, and hard band 
Chandra data. However, we could not identify statisti¬ 
cally significant diffuse X-ray emission in either radial 
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Fig. 3. — Bolometric X-ray luminosity in the (0.05 — 0.15)r200 (left panels) and in the (0.15 — 0.30)r2oo (right panels) region as a function 
of the stellar mass. The upper panels show the 10 — 90 percentile range of the predicted X-ray luminosities in the Illustris Simulation for all 
simulated galaxies within the stellar mass range of (2 x 10^^ — 10^^ Mq. For each stellar mass bin the median X-ray luminosity is marked 
with the horizontal bars. The small points in the middle panels depict the individual predicted X-ray luminosities of all simulated galaxies. 
In the lower panel, triangles show the simulated X-ray luminosity of 42 massive spiral ga laxies from Illustris. The large filled circles points 
show the observed X-ray luminosities for NGC 1961 and NGC 6753 JBogdan et al.ll2013alh along with the 3c7 upper limits presented in this 
work. The large filled diamonds with horizontal bars represent the 3cr upper limits on the X-ray luminosity derived from the stacked data 
of a sub-sampl e of the six most near by galaxies. The uncertainty in the observed X-ray luminosity for NGG 6753 is consistent with the size 
of the symbol JBogdan et al.ll2013al') . Note that there is a large scatter in the predicted X-ray luminosity at stellar masses > 10^^ Mq. The 
scatter and the nearly constant predicted X-ray luminosity for galaxies with stellar mass above 10^^ Mq is the consequence of powerful 
radio-mode AGN feedback, which ejects a significant amount of hot gas from the halos. 
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20^8). The uncertainty in the observed gas mass for NGC 1961 and NGC 6753 is consistent with the size of the symbols 


2013al) . The large filled diamonds with horizontal bars represent the 3a upper limits on the X-ray luminosity derived from 


the stacked data of a sub-sample of the six most nearby galaxies. 


bin. As before, we increased the signal-to-noise ratios by 
stacking the X-ray photons in the three energy ranges. 
The stacked data also did not reveal a statistically sig¬ 
nificant signal. Interestingly, for NGC 266 we detect a 
moderately significant (~ 2cr) X-ray emitting component 
in the (0.05 — 0.15)r2oo region, which possibly originates 
from its hot gaseous corona. We detect 137 ± 60 net 
counts in the 0.3 — 2 keV energy range, where the er¬ 
ror range refers to the statistical uncertainties. Given 
the low signal-to-noise ratio, the importance of system¬ 
atic uncertainties associated with the background sub¬ 
traction is non-negligible. Specifically, depending on the 
background subtraction method, the obtained signal ex¬ 
hibits somewhat different l.Scr — 2.5tT) statistical sig¬ 
nificance. Despite the low significance detection of the 
diffuse emission around NGC 266, we note that the ob¬ 
served count rates a re broadly cons i stent w ith our earlier 
results presented in iBoedan et al.l (|2013bD - for further 
discussion see Section |6l 

These evidences, both individually and combined. 


point out that statistically significant diffuse X-ray emis¬ 
sion is not detected beyond the optical radii (> 20 kpc) 
of the sample galaxies. 

4.2. Constraining the gas parameters 

In the absence of statistically significant detections of 
the extended emission beyond the optical radii, we com¬ 
pute 3cr limits on the characteristics of the hot X-ray 
coronae. To this end, we derive the 3cr upper limits on 
the X-ray counts in circular annuli with r adii of (0.05 — 
0.15)r2 nn and (0.15 — 0.30)r2oo following iBoedan et al.l 
(l2013all . We utilize all available Chandra detectors to 
maximize the number of detected counts. To convert 
the counts to physical units, we assume an optically-thin 
thermal plasma emission (apec) model with kT = 0.2 
keV temperature, 0.1 Solar abundance, and Galactic col¬ 
umn density. To account for the vignetting effects and 
the different sensitivity of various detectors, we employ 
the tailor-made response files for each observation. The 
upper limits on the X-ray counts were converted to the 
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Fig. 5. — The left and right panels shows the temperature and metal abundances of the hot gas within the (0.05 — 0.15)r2oo region as 
a function of the stellar mass, respectively. Both the gas temperature and the abundance refers to luminosity-weighted quantities. The 
vertical bars show the 10 — 90 percentile range, and the horizontal bars show the median values for each stellar mass bin. The observed 
values for NGC 1961 and NGC 6753 are shown as a comparison. Note that to date only these two galaxies offer a measurement of the gas 
temperature and abundance beyond their optical body. 


normalization of the APEC model using: 

^ J ’ 

where D is the distance of the galaxy, while rig and uh 
are the electron and hydrogen number densities, respec¬ 
tively. Based on the upper limits on the emission mea¬ 
sure (f neTiffclV), we compute upper limits on the X-ray 
luminosity, average gas density, and confined gas mass 
within the volume that corresponds to the radial ranges 
of (0.05 — 0.15)r2oo and (0.15 — 0.30)r2oo- Note that 
in these calculations we assumed that the hot gas has 
a spherically symmetric distribution and constant gas 
density within the given bins. Additionally, we derive 
lower limits on the cooling time, which is computed as 
Aooi = (3fcT)/(neA(T)), where A(T) is the cooling func¬ 
tion. Finally, we also derive upper limits on the cooling 
rates of the gas as M = Mgas/tcooh which allows us to 
probe whether cooling of the hot gaseous coronae can 
balance the ongoing star-formation in the disk (Section 

El). 

To place more stringent constraints on the properties 
of hot coronae, we stacked the X-ray photons along with 
the corresponding exposure maps for the six most nearby 
galaxies in our sample, that is for NGC 1097, NGC 2841, 
NGC 5005, NGC 5170, NGC 5529, and NGC 5746. The 
stacked data allow us to compute the average Scr lim¬ 
its on the gas parameters for spiral galaxies in the stel¬ 
lar mass range of (0.7 — 1.4) x 10^^ Mq. Thanks to 
the increased signal-to-noise ratios, these limits are more 
stringent than those obtained for individual galaxies. To 
derive the 3cr limits on the gas characteristics for the 
stacked data, we utilized the median distance, virial ra¬ 
dius, and column density of the sub-sample, and followed 
the procedure described for individual galaxies. The lim¬ 
its are given in Table E] 

Although the main focus of this work is to probe the 
characteristics of hot coronae around lower mass spi¬ 
rals, we extend the present sample by our earlier detec¬ 
tions of the coronae around NGC 1961 and NGC 6753 


(|Bogdan et ahllMlSall . Including these galaxies allows 
us to confront the Illustris Simulation with observational 
results in the stellar mass range of (0.7 —4.2) x 10^^ Mq. 

5. COMPARISON WITH ILLUSTRIS 
5.1. The Illustris Simulations 

The Illustris Project (http://www.illustris-project.org) 
comprises a series of cosmological hydrodynamic simula¬ 
tion runs in a (106.5 Mpc)^ volume performed w ith the 
moving mesh code AREPO (ISpringel et al.l [2M^ . The 
simulations were carried out at a series of resolutions, of 
which the highest resolution, Illustris-1 (or simply Illus¬ 
tris), has a mass resolution of todm = 6.26 x IO^Mq and 
Wbaryon — 1-26 X IO^Mq for the dark matter and bary- 
onic components, respectively. At z = 0 the co-moving 
plummer equivalent gravitational softening lengths are 
1.4 kpc for the dark matter and 0.7 kpc for baryonic col- 
lisionless particles, while for gas it is adaptively set by 
the cell size but cannot be lower than 0.7 kpc. In this 
paper, we rely on the results of the Illustris run. The 
Illustris Simulation discussed here adopts the cosmolog¬ 
ical parameter s that are consistent with the WMAP-9 
measurements (jHinshaw et al.ll20T^ . 

Illustris incorporates the crucial physical processes 
that are i ndispensable to describe ga l axy formation an d 
evolution (|Vogelsberger et ^l2014a 0; iGenel et al.ll2014l) . 
Specifically, it includes primordial and metal line cool¬ 
ing, describes the stellar evolution and the correspond¬ 
ing feedback processes, traces the chemical enrichment 
processes by modeling nine elements, seeds supermassive 
black holes and follows their evolution through accretion 
and mergers, and models the energetic quasar-mode and 
radio-mode feedback of supermassive black holes. De¬ 
tails of the incorp orated physics modules were discussed 
in pre vious works (IVogelsberger et aT1l2013t iTorrev et al.l 
[Ml . A major consequence of the implemented stellar 
and AGN feedback is that the Illustris simulation can 
produce a stellar mass - halo mass relation, which is 
in good agreement with that observed from abundance 
matching studies. Additionally, Illustris can successfully 
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reproduce the observed cosmic star-formation density, 
and can result in a realistic galaxy population. Over¬ 
all, Illustris offers an ideal framework to confront the 
observed and simulated properties of hot X-ray coronae. 

Following our observational approach, we derive the 
characteristics of the hot gas in Illustris in two radial 
ranges. For the inner and outer regions, we utilize the 
0.05r2oo < r < 0.15r2oo and 0.15r2oo < r < 0.30r2oo 
regions, respectively. Since we aim to study the hot gas, 
we only included gas particles, which have a temperature 
of at least 1.16 x 10® K. To derive the X-ray luminosity 
and other p roperties of the gas w e follow the procedure 
described in iNavarro et al.l (|1995ll . 

5.2. X-ray luminosities and gas masses 

Since the X-ray luminosity and gas mass of hot coronae 
have been widely used to test galaxy formation models, 
we first confront these quantities with those predicted 
by Illustris. To facilitate the comparison with previous 
and subsequent works, we perform the comparison in the 
(0.05—0.15)r2oo and in the (0.15—0.30)r2oo radial ranges. 

When investigating the full sample, Illustris predicts 
large scatter in the X-ray luminosity and gas mass at 
every stellar mass (Figs. O S]). Moreover, the median 
X-ray luminosities and gas masses show an increasing 
trend in the stellar mass range of ( 0.2 — 1 ) x 10 ^^ Mq, 
while in M* = (1 — 6 ) x 10 ^^ Mq these values decline 
or remain broadly constant. The reason for the rela¬ 
tively low luminosities and gas masses at M* > 10 ^^ Mq 
is that the radio-mode AGN feedback bec omes domi¬ 
nant a bove these masses. As discussed in iGenel et al.l 
(j2014l) , the radio-mode feedback was tuned in a way that 
it can effectively suppress the star formation in mas¬ 
sive halos, which in turn blows out a significant frac¬ 
tion of the gas from these massive halos, thereby re¬ 
sulting in less massive and less luminous X-ray coro¬ 
nae. For galaxies with less massive black holes, such 
as those residing in low-mass ellipticals or spirals, the 
radio-mode feedback is not dominant due to the weaker 
outbursts. This suggests that spiral galaxies in Illus¬ 
tris host more luminous X-ray coronae than ellipticals. 
To directly probe this possibility, we compare the me¬ 
dian simulated X-ray luminosities for the “textbook” spi¬ 
ral and elliptical galaxies in the stellar mass range of 
(0.7— 1.4) X 10^^ M 0 . Since the median stellar and total 
halo mass of both sub-samples are identical, the different 
X-ray luminosities cannot be attributed to differences in 
the depth of the gravitational potential well. We find 
that within the (0.05 — 0.15)r2oo region the median lumi¬ 
nosity for “textbook” spirals (Lboi = 4.8 x 10^® erg s“^) 
exceeds the median luminosity of “textbook” ellipticals 
(Tboi = 1-1 X 10^® erg s“^) by factor of ^ 44. This 
conclusion contradicts observational results. Indeed, due 
to their more luminous nature, gaseous c oronae around 
mass i ve ellipticals are routinely observed (lForman^ t_alJ 
1985t[O^Sulliyan et al.ll200lt iMathews fc Brighentill2003t 

Bogdan fc Gilfanovll2011ll . as opposed to the significantly 

fainter and less explored coronae around spirals. 

In the (0.05 — 0.15)r2oo region, the observed 3cr upper 
limits on the X-ray luminosities and gas masses are at 
the upper end of the simulated values when considering 
the full Illustris galaxy population. When focusing on 
the 42 simulated “textbook” spirals, we conclude that 


the observed 3cr limits are at the upper end of the model 
predictions. Specifically, the observed limits on the lumi¬ 
nosities and gas masses for the stacked sub-sample exceed 
the median simulated values for the “textbook” spirals. 
In the (0.15 — 0.30)r2oo range the median predicted lu¬ 
minosities and gas masses fall significantly short of the 
predicted upper limits. Based on these, we conclude that 
the Illustris predictions do not contradict the observa¬ 
tional results for spiral galaxies in the stellar mass range 
of (0.7 — 2.0) X 10^^ Mq. However, the predicted lu¬ 
minosities and gas masses are significantly lower than 
those observed for NGC 1961 and NGG 6753, which can 
be attributed to the overly powerful radio-mode feedback 
present in Illustris. 

5.3. Gas temperature and abundance 

Although the non-detections of hot coronae do not 
allow us to constrain the gas temperatures and metal 
abundances, we utilize our previous meas urements for 
NGG 1961 and NGG 6753 ilBogdau et all To 

this end, we confront the predicted luminosity-weighted 
temperatures and abundances with the observed values 
within the (0.05 — 0.15)r2oo region. Since the 42 sim¬ 
ulated “textbook” spirals are less massive than the ob¬ 
served galaxies, we focus on the entire Illustris popula¬ 
tion. 

While galaxies with M* < 10^^ Mq are expected to 
have 0.1 — 0.2 keV gas temperatures, the predicted tem¬ 
peratures shows a rapid increase for more massive galax¬ 
ies. The observed gas temperatures for NGC 1961 and 
NGC 6753 are kT ^ 0.6 — 0.7 keV, which agrees with 
those predicted by the Illustris simulation. 

The right panel of Fig. [5] reveals that Illustris pre¬ 
dicts ~ 0.01 — 0.05 Solar median metal abundances for 
all galaxies in the studied stellar mass range. While the 
observed values exceed the simulated ones by a factor 
of a few, we do not consider this difference robust due 
to the notable systematic unc ertainties associated with 
the abundance measurements (jBogdan et al.l[20T3^ . To 
further constrain the major metal enrichment processes, 
deep X-ray observations are required, which can provide 
the sensitivity to account for the systematic uncertain¬ 
ties in the abundance measurements, thereby constrain¬ 
ing the metal enrichment processes in greater detail. 

5.4. Electron density profiles 

To further probe Illustris, we confront the pre¬ 
dicted electron density profi les with those observ ed for 
NGG 1961 and NGG 6753 (IBogdan et alJ[2nT3^ . The 
right panel of Fig. | 6 ] depicts 61 Illustris galaxies with stel¬ 
lar and virial masses in the range of(3.0 — 4.5)xl0^^ Mq 
and (0.5 — 4.2) x 10^® M 0 , respectively. These values are 
in agreement with those of NGC 1961 and NGC 6753. In 
the right panel of Fig. [51 we compare the observed den¬ 
sity profiles with the predicted density profiles of the 42 
simulated “textbook” spiral galaxies. These galaxies, on 
average, are less massive than NGC 1961 and NGC 6753. 

When comparing the density profiles of the 61 mas¬ 
sive simulated galaxies with those of NGC 1961 and 
NGC 6753, we conclude that the normalization and 
shapes of the profiles are broadly similar. However, 
the agreement is not perfect, because the average nor¬ 
malization of the simulated profiles are somewhat lower 
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Fig. 6.— Electron density profiles as a function of radius. In the left panel, the light grey lines show the simulated profiles for the 61 
galaxies in the Illustris simulation, whose stellar masses is in the range of (3.0 — 4.5) X 10^^ Mq. In the right panel, the light grey lines 
represent the 42 massive spiral galaxies. Note that the stellar mass of these textbook spirals is in the range of (0.3 — 3.2) X 10^^ Mq, and 
hence remain below the stellar masses of NGC 1961 or NGC 6753. The thick red and blue lines represent the observed density profiles 
for NGC 1961 and NGC 6753, respectively. The sho rt dashed red and blue lines show the extrapolation of the density profiles based on 
the best-fit modified /?-models (IBogdan et al.l[2013bli . The regions < 0.02r200 are not plotted, as they are not well resolved close to the 
gravitational softening scale, and they are not sampled by a large number of particles, resulting in larger uncertainties. The radius is 
measured in units of r200. 


than observed for NGC 1961 and NGC 6753, and beyond 
^ 0 . 1 r 2 oo the slope of the simulated profiles are steeper 
than the observed ones. As a caveat, we emphasize that 
beyond ^ 0.15r2oo the density profiles for NGC 1961 and 
NGC 6753 are extrapolated, and hence carry significant 
uncertainties. Therefore, it is not clear how much weight 
should be given to the fact that the observed and simu¬ 
lated profiles exhibit somewhat different slopes at large 
radii. 

A similar picture is obtained when we probe the elec¬ 
tron density profiles of the 42 simulated “textbook” spi¬ 
rals. Interestingly, these profiles have similar normaliza¬ 
tions to those of NGC 1961 and NGC 6753, despite the 
fact that the Illustris galaxies have lower average stellar 
mass. Beyond ~ 0.1r2oo the simulated profiles decline 
more steeply, which may be the combined effect of the 
lower gas mass in these halos and the powerful radio¬ 
mode AGN feedback, which ejects the hot gas to larger 
radii and/or removes it from the dark matter halos of 
galaxies. 

6. DISGUSSION 

Since we are studying spiral galaxies with star- 
formation rates in the ^1 — 6 Mq yr“^ range, it is 
important to probe if these galaxies could be prone to 
starburst driven winds at the present epoch. Indeed, 
theoretical calculations indicate and observational re¬ 
sults confirm that in galaxies with the surface area spe¬ 
cific supernova rate of > 40 SN Myr“^ kpc~^ the detec- 
tion of extraplanar emission is expected (| Strickland et all 
120041) . Note that the surface area is defined as A = 
where D 25 is the diameter of the isophote having the sur¬ 
face brightness of 25 mag arcsec“^. To derive the spe¬ 
cific supernova rates, we compute the supernova rate as 
Rsn = 0.2Lftr/ 10^ ^ Lf7^, where Tftr is the total far- 
infrared luminosity (|Heckman et al.l[l990D . and use the 
D 25 radii given in Table [1] We find that the supernova 
rates are in the range of 0.005 — 0.068 yr“^, and the ob¬ 


served area specific supernova rates are in the range of 
^2 — 36 Myr“^ kpc“^ for all galaxies but NGC 5005 
and ESO 445-081. This implies that we do not expect 
extraplanar gaseous emission driven by the energy input 
of supernovae for the six galaxies in our sample. While 
the specific supernova rates for NGC 5005 and ESO 445- 
081 exceed the critical rate, they remain below the val¬ 
ues o btained for well-known starburst galaxies, such as 
M82 (|Strickland et al.l 120041) . Therefore, while we may 
expect an outflow in NGC 5005 and ESO 445-081, these 
are likely less prominent than those in starburst galaxies, 
hinting that the hot gaseous coronae are not dominated 
by hot gas originating from starburst driven winds. As a 
caveat, we mention that at earlier epochs these galaxies 
likely underwent more active evolutionary phases char¬ 
acterized by gas-rich (minor) mergers, which possibly 
brought a significant supply of cold gas, thereby trig¬ 
gering starburst driven winds. These starburst driven 
winds at high redshifts presumably played a role in lift¬ 
ing a certain amount of hot gas to larger radii, as well as 
enriching the hot gas originating from infall. 

From the derived upper limits on the gas mass and 
the lower limits on the cooling time, we obtain upper 
limits on the cooling rates as M = A/gas/Aooi- Using 
this estimate we probe if the cooling of the hot corona 
can balance the ongoing star-formation in the disk. We 
find that the combined upper limits on the cooling rate 
in the (0.05 — 0.30)r2oo exceed or are comparable to the 
observed star-formation rates for all galaxies in our sam¬ 
ple. Additionally, the cooling of the hot gas in the inner¬ 
most regions (< 0.05r2oo) may also add a contribution, 
since this gas is expected to be denser, and hence has 
a shorter cooling time. This implies that cooling of the 
gas in the corona may be able to provide the cold gas re¬ 
quired for the ongoing star-formation. We note that this 
result is in contrast with that observed for NGG 1961 
and NGG 6753, which systems host significantly hotter 
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gas with cooling times that exceed the Hubble-time. 

Although the Chandra observations of NGC 266 did 
not reveal a statistically significant detection of a hot 
gaseous corona, this result is not in conflict with our 
earlier detection based on ROSAT observations. Specif¬ 
ically, within the (0.05 — 0.15)r2oo region, the observed 
ROSA T count rate was ^ 3.4 x 10“^ s“^ (j Bogdan et al.l 
l2013b[l . As opposed to this, the count rate detected by 
Chandra is ~ 1.1 x 10“^ s“^. Taking into account the ef¬ 
fective area of the ROSAT PSPC and the Chandra ACIS- 
I detectors, we conclude that the factor of ~ 3 difference 
in the count rates may be due to the low temperature 
of the hot gas. Specifically, our observational results are 
consistent with a picture, in which the hot gas corona 
surrounding NGC 266 has a temperature of kT ^ 0.15 
keV. Thus, the non-detection of the hot corona around 
NGC 266 in the deep Chandra observations is likely due 
to the combined effects of the low ACIS-I effective area 
below 0.5 keV energy and the low gas temperature. 

Since in this work - in the absence of statistically signif¬ 
icant detections - we assumed a spectral model to derive 
upper/lower limits on various parameters of the hot gas, 
we briefly overview the uncertainties that may be asso¬ 
ciated with the chosen spectral model. A notable source 
of uncertainty arises from the unknown gas tempera¬ 
tures. For example, if the gas temperature is factor of 
two lower/higher than the assumed 0.2 keV, that would 
imply a factor of about two higher/lower counts-to-flux 
conversion for the ACIS-S detector and factor of about 
three difference for galaxies observed with the ACIS-I 
detector. Additionally, the unknown metal abundances 
may also significantly influence the conversion. For ex¬ 
ample, if the metal abundance are 5 times higher/lower 
than the assumed 0.1 Solar abundance, that would yield 
a factor of ~ 2 increase/decrease in the gas densities and 
gas masses. 

Since in the present study we could not identify a 
so-far unexplored hot corona, the number of galax¬ 
ies with detected coronae remains fo ur: NGC 1961, 
UGC12591, NGC 6753, and NGC 266 (lAndersouTOl 
120111: iDai et all 120111 : iBogdan et al.ll2013a] lbfl. The dif- 
ficulty in identifying hot coronae around spiral galaxies 
with present-day telescopes has two major sources. First, 
the detectable X-ray surface brightness of the coronae is 
low due to their relatively faint but extended nature. 
As a result, the observed signal-to-noise ratios, particu¬ 
larly for lower mass galaxies, do not allow to unambigu¬ 
ously detect these coronae. Second, as pointed out by 
Illustris, the gas temperatures rapidly drops for galaxies 
with M* ^ 3 X 10^^ Mq. However, Chandra and XMM- 
Newton have relatively low effective area at energies be¬ 
low ^ 0.5 keV, and hence are not sufficiently sensitive 
to thermal emission with ~ 0.1 — 0.2 keV. Additionally, 
for such low temperatures the hot gaseous emission is 
prone to absorption, which in turn can drastically re¬ 
duce the observed X-ray flux. The combination of these 
effects makes it difficult to detect hot coronae around 
lower mass spiral galaxies, particularly at radii exceed¬ 
ing ~ 0.15r2oo- 


7. CONCLUSIONS 

In this paper, we studied the hot gaseous coronae 
around 8 normal spiral galaxies, whose stellar mass is 
in the range of (0.7 — 2.0) x 10^^ Mq, using Chandra 
X-ray observations. To confront the observations with 
galaxy formation models, we utilized the hydrodynam- 
ical cosmological simulation, Illustris. Our results are 
summarized below: 

1. We did not detect a statistically significant hot 
corona around any of our sample galaxies in the 
(0.05 — 0.15)r2oo or (0.15 — 0.30)r2oo radial ranges. 
However, for NGC 266 we identified the presence 
of a moderately significant (^ 2a) extended X-ray 
emitting component in the (0.05 — 0.15)r2oo region, 
which emission may originate from hot ionized gas. 
This results is in agreement with our earlier detec - 
tion based on ROSAT data (iBog-dan et al.l[^13bll . 

2. In the absence of statistically significant detections, 
we derived 3 ct limits on the X-ray luminosity, the 
confined gas mass, the electron densities, and the 
cooling times. To derive these limits, we assumed 
that the hot gas is spherically symmetric, has a 
temperature of kT = 0.2 keV and metal abun¬ 
dances of 0.1 Solar. 

3. The derived upper/lower limits, complemented 
with the gas parameters measured for NGC 1961 
and NGC 6753, were confronted with results of 
the Illustris Simulation. We concluded that the 
observed 3cr upper limits on the X-ray luminos¬ 
ity and hot gas mass exceed or are at the upper 
end of the model predictions. For NGC 1961 and 
NGC 6753, we found that the gas temperatures, 
metal abundances, and electron density profiles 
broadly agree with the model predictions. These 
results imply that the observed properties of X-ray 
coronae around spirals are broadly consistent with 
those predicted by Illustris. However, for galaxies 
with massive black holes (mostly giant ellipticals) 
Illustris predicts overly powerful radio-mode feed¬ 
back, which in turn results in under luminous X-ray 
coronae in these systems. 
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